Spin susceptibilities, spin densities and their connection to spin-currents 
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We calculate the frequency dependent spin susceptibilities for a two-dimensional electron gas 
with both Rashba and Dresselhaus spin-orbit interaction. The resonances of the susceptibilities 
depend on the relative values of the Rashba and Dresselhaus spin-orbit constants, which could be 
manipulated by gate voltages. We derive exact continuity equations, with source terms, for the spin 
density and use those to connect the spin current to the spin density. In the free electron model 
the susceptibilities play a central role in the spin dynamics since both the spin density and the spin 
current are proportional to them. 

PACS numbers: 73.63.-b,85.75.-d 



CO 



c3 



i 

-a 
c 
o 
o 



> 

m 
in 

\o 
o 
^t- 

o 

-I— > 



-a 

c 
o 
o 



% 



I. INTRODUCTION 

The ability to manipulate spin states in semiconduct- 
ing, and metallic, nanostructures is now the focus of 
much experimental and theoretical attention. Some spin- 
tronics application are already in usei* 2 " and more have 
been proposed, ranging from spin FET's 3 ^^ to spin 
qubits^i In many cases the spin-orbit interaction can 
be used to manipulate the electron spin via electronic 
means £ This sort of electronic manipulation is impor- 
tant since technologically, electric field control of spins 
is preferred over magnetic held control. Also, electronic 
control of spins has revealed many interesting physics in 
experiment involving semiconductor heterostructuresAiS 

Already quite some time ago, it was proposed that 
electric fields could lead to magnetization in anti- 
ferromagnetic materials fulfilling certain symmetries^ 1 ", 
this is the so called magnetoelectric effect. The 
magnetoelectric effect was considered for conductors 
with special symmetry properties^ and for pyro- 
electric superconductors^ 3 -. Electric field induced spin- 
orientation in semiconductors due to linear in momen- 
tum spin-orbit interaction was also discussed in Refs. 
Il4lll5lll6l There has been renewed interest in this subject 
since this induced spin polarization might serve as spin- 
injectors in certain semiconductor heterostructuresiiiiS 
Only recently there have been experiments which seem 
to demonstrate such current induced polarization of spins 
via spin-orbit interaction^*^ 

A somewhat related effect is the so called spin-Hall 
effect. In the normal Hall effect the electrons are de- 
flected by the Lorentz force, but in the spin-Hall ef- 
fect spins are scattered by impurities preferentially to 
the left (right) if their spins point 'up' ('down') due to 
spin-orbit interaction. 21 ! 22 This will result in a spin cur- 
rent, but no net charge current, flowing perpendicular to 
the applied charge current. Such an extrinsic transverse 
spin current was already investigated by Dyakononov and 
Perel some time ago in Ref. |23| and more recently by 
Hirsch 24 . Presently there is much interest in spin-orbit 
mediated spin-Hall effect in semiconducting hcterostruc- 
tures. Here the effect can appear in holo 25 i 26 i 27 ' 28 i 29 i 30 i 31 
or electron 32 i 33 i 34 i 35 i 36 i 37 i 3s i 39 i 40 i 41 i 42 doped semiconduc- 



tors due to band structure properties and impurities are 
not necessary. Thus, the term intrinsic spin Hall ef- 
fect is used to distinguish this from the previously dis- 
cussed mechanism since it occurs already in the absence 
of impuritiesi 2 ^ 2 " Due to differences in the band struc- 
ture of the holes and electron there is an important dis- 
tinction between the two since in the electron system the 
spin Hall conductance (the ratio of the transverse spin 
current and the applied electric field) reaches a universal 
value of e/87r in clean systems^ However, impurities are 
believed to modify this universal valu o 33 i 40 and in addi- 
tion, the exact influence of vertex corrections on the spin- 
Hall conductivity is currently under investigation 3 ^ 4 ^ 3 .. 
In this paper we consider the spin susceptibilities of 
a 2DEG with Rashba and Dresselhaus spin-orbit cou- 
pling. Due to the strong 2DEG confinement the Dres- 
selhaus coupling reduces to terms linear in momentum. 
For such linear momentum spin-orbit coupling the spin 
susceptibility can be used to characterize other transport 
properties. We calculate the susceptibilities using a free 
electron model and relate them to the electric field in- 
duced spin density. Also, we derive continuity equations 
(with source terms) for the spin density and spin cur- 
rent, similar to the equations already derived for only 
Rashba coupling^ These equations are exact operator 
identities and via them we can relate the spin current 
to the spin density. Via these relations the spin current 
(which is non-trivial to measure) can be connected to the 
spin density, or magnetization, which is easier to detect. 
The susceptibilites play a central role in the free electron 
model, since the Fourier transform of the spin density and 
the spin current are proportional to the susceptibilities. 



II. THE MODEL 

We consider a two dimensional electron gas (2DEG) 
with Rashba and Dresselhaus spin-orbit coupling. In the 
absence of external fields the Hamiltonian may be written 

as 
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vl 
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-(PyCTx -Px°y) + TiPyVy 



Px°z 



(1) 



where a and (3 are the Rashba and Dresselhaus (linear) 
coefficient in a 2DEG. It is easy to see that [p x ,H] = 
[p y , H] = 0, and thus we seek eigenstates of the form 



ipk,s( r ) — (r\ks) 



ik-r 



u s (k) 



(2) 



where U s (k) is spinor to be determined and A is the sys- 
tem area. Since the momenta are conserved it is possible 
to introduce an effective magnetic field 



The eigenf unctions in Eq. Q have the interesting prop- 
erty that they depend on the spin-orbit coupling parame- 
ters a and (3 only via the angle (p. This allows one to take 
the limit a, (3 — > such that the angle <f> remains fixed and 
the resulting eigenvectors in Eq. (JJJ are also (degenerate) 
eigenvectors of the free electron Hamiltonian. 44 As was 
pointed out in Ref. y3, the Kramers conjugate state of 
i>k,s(i') is ip-k tS (r), i.e. they belong to the same branch. 
III. SPIN SUSCEPTIBILITIES AND CURRENT 
INDUCED MAGNETIZATION 



-(3k x + ak y 

(3k y - ak x 





(3) 



The eigenspectrum can be written in terms of this effec- 
tive magnetic field, the eigenenergies of the Hamiltonian 
being (s = ±1) 






and the corresponding spinor 



(4) 



For a weak driving field the response of the system is 
obtained by the Kubo formalism. Due to the spin-orbit 
coupling, a pure electric field driving results in a non-zero 
magnetic response. Since the spin-orbit term in Eq. Q) 
is linear in momenta, both the response functions due 
to magnetic and electric perturbation can be expressed 
with spin susceptibilities. The r\ = x,y, z component of 
the spin density operator is defined as 



pvM =y^ n 
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S(r - r„), 



(9) 
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where T + = T x +iT y . Note that Eqs. (gj and J5J are also 
valid for an in-plane magnetic field B using the substitu- 
tion r — ► r + ig/i^JB, where g and [Ib are the effective 
^-factor and Bohr magneton, respectively. 

Writing the quasi-momentum in polar coordinates fc = 
/c(cos#,sin#), the eigenspectrum of the Hamiltonian in 
Eq. QJ becomes 



E s (k) 
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+ S ky/( a 2+f32)g(6), 



(6) 



where g{0) = 1 — sin(2</>) sm(26) determines the 
anisotropy of the Fermi surfaces and the corresponding 
eigenspinors are 



1 



u s (k) = — 



(cos(0)e l -\-isin((f>)e z ) 



(7) 



Here we have introduced the following parameterization 
of the spin-orbit coupling strength 



sin((/>) 



V^+^ 
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where r„ and a n ^ are the position operator and Pauli 
matrix, respectively, of the nth electron. For a transla- 
tionally invariant system, the wave vector and frequency 
(q, uj) dependent susceptibilities are 



W(g,«)=/ dte^—dpniq^^n'i-q)}}, (10) 

(11) 
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where fk s — f(E s (k)) with / being the Fermi distri- 
bution function and we used the notation (. ..)ha = 
(ks\ . . . \ks). The frequency should be viewed as iu> — » 
i(uj + ih~ l e) to regularize the integral. The operators in 
Eq. I|lll) refer to single particle operators. This suscepti- 
bility is a spin density response function and to get the 
magnetization response function, each spin density oper- 
ator in Eq. I|llf) should be multiplied with the electron 
effective magnetic moment 5/13/2. Using Eq. (|11|) and 
the eigenspectrum represented by Eqs. JBJ) and J7J) the 
susceptibilities for a spatially homogeneous perturbation 
(q = 0) become 
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Note that all q = dependence has been dropped for 
clarity. Here we assumed zero temperature and the Fermi 
distribution function was replaced by a step function. 
The resulting Fermi contours k±(9) are the solutions of 



k% = k±(e) 2 ±k±(e)k S oVg(eji 



(14) 



where kp = 2mep/K 2 is the squared Fermi momentum 
and kso — m^f a 2 + B 2 /h 2 . The k integral results in 
a linear term and a term involving an inverse tangent 
in k±(9) and the subsequent angular integrals cannot be 
solved analytically. The inverse tangent can be expanded 
in powers of kso /kF ^ 1? resulting in the following low- 
est order result 



Xxx(uj) = Xyy(u) 



2irh 2 



( 1+ n; 
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(15) 

where the resonance energies are e\ — 8eso£f(1 ± 
sin(2</>)), with £50 = m(a 2 + (3 2 )/h 2 . Using the same 
procedure we can calculate the off-diagonal susceptibili- 
ties in a similar manner 



1 \ ( \ 5 x{u)~Xxx{u) . .,,,-, , 
Xx V \u) = Xyx(u) = „ : „ /ojA sm(2<p)dx(u) 
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where we have defined 

n 

<*x( w ) = 
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(16) 



(17) 



The magnetization is related to the spin density in Eq. @ 
through m(r) = ■kgPBPi'r) which leads to the standard 
linear response relation 

m v (q,u) = y-gp B ) x w '(g, u)B v ,(q, u). (18) 

To obtain the Pauli paramagnetic susceptibility one 
should take the following order of limits^ 5 * 



lim lim XvV («?> w ) 

q—>0 u^O 



m 



(19) 



The diamagnetic contribution can be disregarded since 
we assume an in-plane magnetic field. The susceptibili- 
ties in Eqs. (|12f> and l|13l) are calculated for the reverse 
order of limits done in Eq. (|19fl . These lj ^ sus- 
ceptibilities are the spin-orbit contribution coming from 
the region in fc-space where only one (s=— 1) of the two 
branches is occupied. 55 Also, the spin current (see discus- 
sion below) which results from the spin-orbit interaction 
is non-zero due to contributions from the fc-space area be- 
tween the two Fermi contours ,22iiS2i Thus we only focus 
on this contribution when we relate the susceptibilities 
to the spin-orbit mediated spin densities and the spin 
currents. 

The value of the resonance frequency is determined 
by y/8£ S o£F ~ 0.16 meV=40 GHz for typical GaAs 
parameters^: a = 0.5 x 10 -9 meV m and electron density 



n e — 4x 10 15 m -2 . For lower frequencies the susceptibil- 
ities remain nominally constant. In the limit a <C B the 
lowest order contribution to Eqs. (fl 51) and (|16fl become 
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Here we have not included impurities and thus the regu- 
larization parameter e can strictly only be attributed to 
an adiabatic turning on of the external electric or mag- 
netic field. 

In the absence of electric and magnetic fields, the spin- 
orbit interaction does not give rise to a net magneti- 
zation. Even though the spin-orbit interaction has the 
form of a momentum-dependent magnetic field, the to- 
tal contribution averages to zero. 47-48 However, although 
the spin-orbit induced splitting does not give rise to an 
equilibrium magnetization, there is an asymmetry in the 
local magnetic field in momentum space, i.e. the local 
magnetic field is odd under wave vector reversal, and any 
translation of the Fermi sphere away from the T-point 
will induce a magnetization^ 5 - Applying a homogeneous 
electric field E(r,t) — E°e~ luJt to the system will give 
rise to the following time dependent perturbation 



V(t) 



_ E o e -iut 



j(q = 0) 



(22) 



where j(q) is the Fourier transform of the current density 
operator 

j(r) = £jn(r) = ]Ti{<5(r-r„),M, (23) 



and the velocity operator v n for the Hamiltonian in Eq. 
||TJ is given by 



rn Pn,x ft 0"' 



S"n>J 






C x 



(24) 



From linear response theory the Fourier transform of the 
electric field induced spin density may be written as 



/ / w e£j x {ui) (a p \ 

{Px{u)) = — ; (tXxj/(w) + -XxxM 

ilu \n h / 

— J^\k xM + k xM ) (25) 

/ 1 w &E x (w) fa (3 \ 

(Pj/M) = t- kfe^J + -rXxyiu)) 

eE v (uj) (a , . 8 , ,\ ,„„. 
it~^\H Xxy{u) + Ji Xxx{uj) )- (26) 

This result is reminiscent of the pure magnetic field in- 
duced spin density, except here eE v {ui)/iu} plays the 
role of magnetic field, via the spin-orbit coupling. As 
was pointed out in Refs. Il5ll4ll the dc limit corresponds 



to replacing the frequency with momentum scattering 
iuj — » — 1/r. 

Multiplying Eqs. (|25|l and H26[l by the sample area 
will give the total number of induced magnetic moments, 
measured in units of gps /2. Applying an electric field 
E sw lOOV/cm to a GaAs 2DEGs with high mobility 
(H/t ~ 10~ 2 meV) and a sample area of A = (500 /im) 2 , 
the number of magnetic moments (Bohr magnetons) 
would be around 2.5 x 10 7 . For a 2DEG thickness of 
a few nm these magnetic moments produce a magnetic 
field of the order 10~ 6 T. Probing 2DEG properties us- 
ing ESR techniques has been succesfully used to de- 
termine the spin-orbit splitting 49 and other 2DEG spin 



propertie; 



,50,51 



Using similar ESR techniques, the spin- 



orbit coefficients a, f3 could in principle be determined 
by measuring different spin density component for dif- 
ferent direction of driving current, as a function of a 
which could be tuned by gate voltages. 46 ' 52 '^ Such an 
induced spin density could be detected by Faraday rota- 
tion measurements jiS 



IV. CONNECTION TO SPIN CURRENT 

The usual way of deriving the operator version of the 
particle continuity equation is to start from the definition 
of the density operator 



p(r,t) = ^<5(r- 



.(*)), 



(27) 



and from there one can derive the standard continuity 
equation 






V.j(r,t)=0, 



(28) 



from the Heisenberg equation of motion for p(r,t). The 
form of the current density in Eq. (|28|l is uniquely de- 
termined by the Heisenberg equation of motion. For the 
Hamiltonian in Eq. Q the current density is given by 
Eq. (|23l) . using the velocity operator in Eq. (|24|) . In the 
case of the spin density operator 



Pn 



( r > *) = X] a ri,ri(t)6{r - r„(t)), 



(29) 



the same procedure of evaluating the Heisenberg equa- 
tion motion will not result in a unique definition of the 
associated spin current density. The reason is that the 
precession due to the momentum dependent spin-orbit 
magnetic field introduces additional terms into the equa- 
tion of motion. In order to proceed one has to postulate a 
form for the spin current density. The most widely used 
definition, and physically reasonable, is the following 



f{r) 



1 



Wn,r,Jn(r)}. 



(30) 



This form of the spin current is Hermitian and reduces to 
the standard spin current form when the velocity oper- 
ator is spin independent, 5 Having determined the form 



of the spin current the resulting continuity equations be- 
come 

p x (r,t)+V-j x (r,t) = —fjr,t) + —ff y (r,t) 



dt 
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U:(r,t)+fJr,t)). 



(33) 



Similar equations have already been derived for pure 
Rashba coupling^ The above Eqs. I|31|) - (|33|l are exact 
relation for a systems with spin-orbit coupling linear in 
momentum and including impurities would not change 
their forr m 5657 

Taking the thermal average of Eqs. l|5T|l - l|55j l gives par- 
tial differential equation connecting the spin densities 
and spin currents. Based on these equations one can 
make a few observation on the nature of the spin cur- 
rent, without explicitly solving them. First of all, for a 
homogeneous system in the stationary limit the rhs of 
Eq. (|33l) must vanish for all values of a, (3. This is triv- 
ially satisfied for all (j v ) = 0, but more interestingly also 
when (j*) = ~{jl) and (j%) = -($}< It is easy to show 
that the latter case is true in equilibrium 
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(34) 


(f x ) = 


-<£> - 


1 a ( m\ 2 , 2 


-p 2 ), 


(35) 



which covers the results of Ref.|33 as special cases. In par- 
ticular, the current expectation values vanish at a = ±/3 
due to the additional concerved quantity arising at these 
points^. Furthermore, Eq. (|33|l shows these equilibrium 
currents do not act as source terms for the spin density, 
since the rhs always vanishes. 

Let us now consider a homogeneous system such that 
the divergence terms vanish. By using Eqs. (|31|l and 
(|32|l one can derive the following identity for the Fourier 
transform of the x and y component of the average z- 
polarizcd spin current 



(£M> 
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h 2 iuj(a(p x (uj)} + (3(p y (uj)}) 
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(37) 



These relations establish a connection between the spin 
current components (j£ y ) and the x, y components of 
the spin density, in the frequency domain. This is quite 
useful since the spin current, which is hard to detect, is 
related to a quantity which is easier to measure. Also, 



Eqs. i|36|) and Q37J I is a good starting point for connecting 
spin current and spin density response functions using 
standard Kubo formalism. 

Let us now assume a homogeneous electric field applied 
in the y direction. The spin conductivities are defined as 
the ratio of the spin current and applied electric field 



>ni 



(«) = 



|0») 

Ey(w) ■ 



(38) 



where the factor h/2 in the definition of the spin conduc- 
tivity is due to our definition of the spin current in terms 
Pauli matrices and not the spin operators, i.e. S = hcr/2. 
Using Eqs. (|25[) and H2bj) to relate the spin density to the 
susceptibilities we obtain the following result for the ac 
spin conductivities 
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(39) 
Xx*(«0). (40) 



These equations for the spin conductivities, along with 
Eqs. (|15|l and (I17|l . in the oj — > limit reproduce the 
result in Ref.|34]up to a sign convention for (3. For either 
pure Rashba or Dresselhaus we have the following result 



Vxviu) 
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(41) 
(42) 



the upper (lower) sign refers to a pure Rashba (Dressel- 
haus). Here one can replace e with the impurity scatter- 
ing broadenin g H/ t, which reproduces the zero frequency 
result in Ref. y3[ Taking the limit e — > (r — > oo) 
recovers the universal limit of spin Hall conductance, 



'.y 



(w-0) 



32 



V. CONCLUSION 

We calculated the frequency dependent spin suscepti- 
bilities for a 2DEG with both Rashba and Dresselhaus 
spin-orbit coupling. The suspectibilities have resonance 
peaks whose position depends on the relative magnitude 
of the Rashba and Dresselhaus coefficients. The position 
of the resonance could be tuned via the gate dependence 
of a. We derive a connection between spin density and 
the spin current starting from the Heisenberg equation 
of motion for the spin density. Unlike the particle den- 
sity, the resulting continuity equations have spin current 
source terms due to the spin-orbit interaction. Using 
these relation we connect the spin current to the spin 
density. In the free electron model the spin current, and 
the spin-Hall coefficient, are proportional to the suscepti- 
bility. Indeed, using the continuity equations and the sus- 
ceptibilities calculated here we recover the universal (bal- 
listic) value of the spin-Hall conductance. Such relations 
(Eqs. (K)(:i[) and (|37|0 might help elucidate the nature of 
spin currents in a similar way to Ref. 41 which discussed 
the relation between the spin current and the dielectric 
function. By calculating the spin density with the correct 
impurity contribution would automatically give the spin 
current. We are confident that these and similar consid- 
erations will contribute to a deeper understanding of the 
role of impurities in the spin-Hall effect. 
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Assuming that VFq <C fcsT <S ^/£so s f, the dominant 
contribution comes from the region between the two Fermi 
contours (the spin-orbit contribution) if huj ^> EF<l/kF- 

56 When Hamiltonians representing non-magnetic impurites 
^imp — ^imp( T ') are added to Eq. Q, the commutator in 
the Heisenberg equation of motion for an operator 0(r, cr) 
becomes [0(r,cr),H + Hi mp ] = [0(r, cr),H], where O is 
any operator that does not explicitly depend on p. Both 
the particle and spin density are such operator and thus 
the equations of motion remain the same in the presence of 
impurities. The complications due to the impurities appear 
in the calculation of the average values of the operators. 

57 Homogeneous (q~0) spin systems driven by a constant 
field E will exhibit damping of magnetization, i.e. 
dt{p v (t — > oo)) = for any finite amount of dissipation. 
Then, from Eq. 1311 and the relation limt^oo dt{px(t)) = 
lmis^o s {px(s)) for the Laplace transform, we get = 
lim 3 _*o s 2 p x (s) — 2m $2 Ea <Jxy(s — > 0), where we used 
(jx{ s )} — ^&xy( s )Eo/s and a — 0. Thus, the dc spin 
conductivity vanishes for any dissipation mechanism. The 
same follows directly from the other magnetization com- 
ponents and in the case of a / 0. 



